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The PSAM (Probabilistic Structural Analysis Methods) program, funded by NASA Lewis 
Research Center, is developing a probabilistic structural risk assessment capability for the SSME 
components. PSAM is currently in the seventh year ofa two-phase, ten-year contract. An advanced 
probabilistic structural analysis software system, NESSUS (Numerical Evaluation of Stochastic 
Structures Under Stress), is being developed as part of the PSAM effort to accurately simulate. 
stochastic structures operating under severe random loading conditions. 

A central part of the NESSUS system is a finite element analysis (FEA) module. FEA is 
generally known to be computer intensive. Thus, the conventional Monte Carlo method, which 
requires a large number of simulations (i.e., a large number of deterministic computer runs), is too 
time-consuming to be practical for probabilistic FEA analysis. One of the major challenges in 
developing the NESSUS system is the development of the probabilistic algorithms that provide 
both efficiency and accuracy. The main probability algorithms developed and implemented in the 
NESSUS system are efficient, but approximate in nature. In the last six years, the algorit hms have 
improved very significantly. 

In probabilistic FEA analysis, a good index for measuring the computational efficiency is the 
number of deterministic solutions required for the user-selected performance function. To minimize 
this number, denoted as M, the first approach taken by PSAM was to generate a response surface 
over a "wide" range (say, ± 3 standard deviations for each random variable). Once the response 
surface is generated, a fast probability integration (FPI) algorithm [Ref. 1] can be used. In practice, 
when the number of random variables is not small, M might be too large, and the FEA pan tends 
to dominate the total computational time. Moreover, the response surface approach does not 
generally provide sufficient accuracy unless an expensive iterative procedure is applied to update 
the response surface at focused regions [Ref. 2]. 

To improve the efficiency, the concept of FPI was applied directly to guide the FEA to develop 
a good approximate performance function. The basic concept is to use the initial information from 
the conventional mean value first order (MVFO) solutions to identify regions that are 
probabilistically more likely for the given performance function values, then move the FEA to these 
regions. MVFO requires (n + 1) deterministic solutions, where n is the number of random variables, 
and provides approximate mean and standard deviation for the performance function. However, 
in probabilistic structural analysis, it is more desirable and often necessary to have knowledge in 
the whole distribution function (CDF). The advanced mean value (AMV) method was developed 
to provide the guidance for the FEA "move," and for efficiently generating performance CDF based 
on the MVFO solution [Refs. 3, 4]. The AMV method has been found to be quite effective for a 
wide variety of engineering problems. Further procedures requiring more M were also developed 
to improve the accuracy of the AMV method. In summary, M - n + 7 is believed to be the minimum 
number required to obtain reasonably accurate probabilistic output that includes the performance 
CDF, and the probabilistic sensitivity factors for the input random variables. 
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Establish CDF Using Fast 
Probability Integration Algorithm 


• Define Limit-State 


> Approximate Performance Function 

At One or More Probability Significant Regions 

> Compute Probability Based on Approximate Functions 


Joint Probability Density 
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ADVANCED MEAN VALUE (AMV) METHOD 

• Conventional Mean Value First-Order (MVFO) method 

First-order Taylor’s series expansion at mean values; 

Z=a.+ZaA(=Z,) 

Valid for small standard deviations. 

• Advanced Mean Value (AMV) First Order Method 

Z* = Z,+//(Z I ) 

Features: 

- //(Z,) introduced to minimize truncation error . 

- Iteration procedure available to find H(Z,). 

- Can be used to detect non-monotonic functions. 
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AMV - Based Iteration Procedure 





124 




AMV EXAMPLE: RANDOM VIBRATION 
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Turbine Blade Verification Example 
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CUMULATIVE PROBABILITY 


MODE 1 - TANGENTIAL FLCX MOOE 
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Fast Convolution Procedure 


1. Dependent Non-normal RVs to Independent Std. Normal RVs. 

2. Find Most Probable Point and Construct Second-order Approx. 

3. Eliminate Product Terms by Orthogonal Transformation. 

4. Transform to Linear Polynomial. 

5. Apply Convolution Theorem. 


Fast Convolution 

• SUM OF INDEPENDENT RANDOM VARIABLES 

Z=X l +X 2 +...+X i + ...+X ll 

• CHARACTERISTIC FUNCTION 

4>((o) = [“ f(x)e iva dx 
d> 2 ((0) = 

. USE FAST FOURIER TRANSFORM TO COMPUTE 4>(oj) 

. USE INVERSE FTT TO COMPUTE PDF OF Z 

• FOR NON-LINEAR Z . USE MPP . QUADRATIC APPROXIMATION AND 
LINEARIZATION 
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Linearization 


ft # ft « • 

« 2 (u) =a 0 + I a, (Mi -U,)+ I Vm, - Mi ) + I I Cy(l4i - U, ) (Uj - Uj ) 
i - 1 f - I I — I > — I 

• QUADRATIC APPROXIMATION AND LINEARIZATION 


^(X)sfl 0 + 

1 f «i 2 

C ° 0 ° 4,?i b t 

g (Y) = c 0 + 

r <=[ x ‘ + itJ =iX ‘ +A ‘ f 


• LOG-TRANSFORMATION 

g(X) s a 0 + £ a, Injr, + ^(Inar,) 2 

i * I 
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Fast Convolution Example 

Z s R - S 






Fast Convolution Validation Problem 


Random variable S has bi-modal pdf: 



Limit state function: 

g = R-S = 0 



Strength: R~Lognormal~(20., 5.) 

Stress: S~Bi-modal~(pi l p2, oi, 02 ) = (10,2,40,2) 


Probability of failure: 


a 

b 

Exact 

Previous Method 

Improved Method 

19 

20 

6.331 E-2 

8.29E-2 

6.285E-2 

99 

100 

2.307E-2 

Numerical Problem 

2.347E-2 


SYSTEM RELIABILITY ANALYSIS METHODOLOGY 



COMPUTE: 

P[(9 1 < 0) U (g 2 < 0) 
U((g 3 <0)n(g 4 <0))] 
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ADAPTIVE IMPORTANCE SAMPLING METHOD FOR 
SYSTEM RELIABILITY ANALYSIS 

METHODOLOGY 



Th# Concept Ol tht Curvature Baaad Adaptlva Importanca Sampling Mtlhod 


Tht Concept of Gant rating Importanca Sampling Potnia 
lor Multlpla umli Stataa 


ADAPTIVE IMPORTANCE SAMPLING METHOD FOR 
SYSTEM RELIABILITY ANALYSIS 


ROTORDYNAMICS EXAMPLE 



Curvatura-BasAd Important# Sampling Point • 

Probability of Inability Sampling Point* lor X, > 0 and >„ > 0 In m# u-Spac* lor X, > u and X, > l> In tft# u-Spac# 

(5000 Mont# Cano Samplaa) 


130 





